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Abstract 

We present an efficient algorithm that computes the Minkowski sum 
of two polygons, which may have holes. The new algorithm is based on 
the convolution approach. Its efficiency stems in part from a property 
for Minkowski sums of polygons with holes, which in fact holds in any 
dimension: Given two polygons with holes, for each input polygon we can 
fill up the holes that are relatively small compared to the other polygon. 
Specifically, we can always fill up all the holes of at least one polygon, 
transforming it into a simple polygon, and still obtain exactly the same 
Minkowski sum. Obliterating holes in the input summands speeds up the 
computation of Minkowski sums. 

We introduce a robust implementation of the new algorithm, which 
follows the Exact Geometric Computation paradigm and thus guarantees 
exact results. We also present an empirical comparison of the performance 
of Minkowski sum construction of various input examples, where we show 
that the implementation of the new algorithm exhibits better performance 
than several other implementations in many cases. In particular, we com¬ 
pared the implementation of the new algorithm, an implementation of the 
standard convolution algorithm, and an implementation of the decomposi¬ 
tion approach using various convex decomposition methods, including two 
new methods that handle polygons with holes—one is based on vertical 
decomposition and the other is based on triangulation. 

The software has been developed as an extension of the 2D Minkowski 
Sums package of Cgal (Computational Geometry Algorithms Library). 
Additional information and supplementary material is available at our 
project page http: //acg. cs . tau. ac. il/projects/rc 
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1 Introduction 


Let P and Q be two point sets in The Minkowski sum of P and Q is 
defined as P ® Q = {p + q\p & P, q £ Q}- In this paper we focus on the com¬ 
putation of Minkowski sums of general polygons in the plane, that is, polygons 
that may have holes. However, some of our results also apply to higher di¬ 
mensions. Minkowski sums are ubiquitous in many fields m, including robot 
motion planning US, assembly planning [7], and computer aided design 0. 

1.1 Terms, Definition, and Related Work 

During the last four decades many algorithms to compute the Minkowski sum 
of polygons or polyhedra were introduced. For exact two-dimensional solutions 
see, e.g., [5]. For approximate solutions see, e.g., m and m- For exact and 
approximate three-dimensional solutions see, e.g., el na, mi, and [21] . 

Computing the Minkowski sum of two convex polygons P and Q is rather 
easy. As P ® Q is a convex polygon bounded by copies of the edges of P and 
Q ordered according to their slope, the Minkowski sum can be computed using 
an operation similar to merging two sorted lists of numbers. If the polygons are 
not convex, it is possible to use one of the two following general approaches: 

Decomposition Algorithms that follow the decomposition approach decom¬ 
pose P and Q into two sets of convex sub-polygons. Then, they compute 
the pairwise sums using the simple procedure described above. Finally, they 
compute the union of the pairwise sums. This approach was first proposed 
by Lozano-Perez [2D]. The performance of this approach heavily depends on 
the method that computes the convex decomposition of the input polygons. 
Flato et al. [T] described an implementation of the first exact and robust ver¬ 
sion of the decomposition approach, which handles degeneracies. They also 
tried different decomposition methods, but none of them handles polygons with 
holes. 

Ghosh [D] introduced slope diagrams —a data structure that was used later 
on by some of us to construct Minkowski sums of bounded convex polyhedra 
in 3D |S]. Hachenberger El constructed Minkowski sums of general polyhedra 
in 3D. Both implementations are based on the Computational Geometry Algo¬ 
rithms Library (Cgal) and follow the Exact Geometric Computation (EGC) 
paradigm. 

Convolution Let Vp = (p 0 , ■ ■ ■ ,p m -i) and Vq = (q 0 ,... ,q n - 1 ) denote the 
vertices of the input polygons P and Q, respectively. Assume that their bound¬ 
aries wind in a counterclockwise order around their interiors. The convolution 
of these two polygons, denoted P*Q, is a collection of line segments of the fonrr] 
\pi+qj,Pi+i+qj\, when the vector PiPi+i lies counterclockwise in between qJZJq$ 

3 Addition of vertex indices is carried out modulo n for P and modulo m for Q. 
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Figure 1: The convolution of a convex and a non-convex polygon; winding 
numbers are indicated in brackets; dotted edges are left out for the reduced 
convolution. 


and qjqj+t and, symmetrically, of segments of the form [p, + q 3 . + qj + i], when 

the vector qJqJ^X lies counterclockwise in between pi_\p\ and piPi + {. 

According to the Convolution Theorem stated in 1983 by Guibas et al. GDJ, 
the convolution P * Q of two polygons P and Q is a superset of the boundary 
of the Minkowski sum P © Q. The segments of the convolution form a number 
of closed (possibly self-intersecting) polygonal curves called convolution cycles. 
The set of points having a nonzero winding number with respect to the convolu¬ 
tion cycles comprise the Minkowski sum P © Q 0 However, this theorem has not 
been completely proven. In the introduction of the thesis of Ramkumar |22| . 
there are statements about the correctness of the Convolution Theorem. Some 
of these statements are still given without proofs. 

Wein [25] implemented the standard convolution algorithms for simple poly¬ 
gons. He computed the winding number for each face in the arrangement in¬ 
duced by the convolution cycles and used it to determine whether the face is part 
of the Minkowski sum or not; see Figure [l] Wein’s implementation is available 
in CGAL [26], and as such, it follows the EGC paradigm. 

Kaul et al. m observed that a segment [p* + qj,Pi +1 + qj] (resp. \p.i + 
qj,Pi + 9j+i]) cannot possibly contribute to the boundary of the Minkowski 
sum if qj (resp. pj) is a reflex vertex (see dotted edges in Figure [TJ . The 
remaining subset of convolution segments, the reduced convolution, is still a 
superset of the Minkowski sum boundary, but the idea of winding numbers can 
not be applied any longer as there are no closed cycles anymore. Instead, Behar 
and Lien [2] , first identify faces in the arrangement of the reduced convolution 
that may represent holes (based on proper orientation of all boundary edges 
of the face). Thereafter, they check whether such a face is indeed a proper 
hole by selecting a point x inside the face and performing a collision detection 
of P and x © — Q. Their implementation exhibits faster running time than 
Wein’s implementation. However, although it uses advanced multi-precision 
arithmetic, it does not handle some degenerate cases correctly. The method 

4 Informally speaking, the winding number of a point p S R 2 with respect to some planar 
curve 7 is an integer number counting how many times does 7 wind in a counterclockwise 
orientation around p. 
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was also extended to three dimensions [T5| . 

Milenkovic and Sacks m defined the Monotonic Convolution, which is an¬ 
other superset of the Minkowski sum boundary. They show that this set defines 
cycles and induces winding numbers, which are positive only in the interior of 
the Minkowski sum. 

1.2 Our Results 

We present an efficient algorithm that computes the Minkowski sum of two 
polygons, which may have holes. The new algorithm is a variant of the algorithm 
proposed by Behar and Lien [2], which computes the reduced convolution set. 
In our new algorithm, the initial set of filters proposed in [2] is enhanced by 
the removal of complete holes in the input. This enhancement reduces the size 
of the reduced convolution set even further. The enhancement is backed up by 
a theorem, the proof of which is also presented; see Section [2] Moreover, we 
show that at least one of the input polygons can always be made simple (before 
applying the convolution). These latter results are applicable to any dimension 
and are independent of the used approach. In addition, roughly speaking, we 
show that every boundary cycle of the Minkowski sum is caused by exactly one 
boundary cycle of each summand; see Section [2j It implies that we can compute 
the convolution of each pair of boundary cycles of the summands separately. 
This result is also applicable to any dimension and it is independent of the used 
approach. However, applying it to the decomposition approach requires the 
ability to handle unbounded polygons. 

We introduce an implementation of the new algorithm. We also introduce 
implementations of two new convex decomposition methods that handle poly¬ 
gons with holes as input—one is based on vertical decomposition and the other 
is based on triangulation. These two methods can be directly applied to com¬ 
pute the Minkowski sum of polygons with holes via decomposition. All our 
implementations are robust and handle degenerate cases. 

We present an empirical comparison of all the implementations above and 
existing implementations; see Section [4] We show that the implementation of 
our new algorithm, which computes the reduced convolution set, exhibits better 
performance than all other implementations in many cases. 


2 Filtering Out Holes 

The fundamental observation of the convolution theorem is that only points on 
the boundary of P and Q can contribute to the boundary of P®Q- Specifically, 
the union of the segments in the convolution P *Q, as a point set, is a super-set 
of the union of the segments of the boundary of P © Q. 

The idea behind the reduced convolution method is to filter out segments 
of P * Q that can not possibly contribute to the boundary of P © Q using a 
local criterion; see Section 0 In this section we introduce a global criterion. 
We show that if a hole in one polygon is relatively small compared to the other 
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polygon, the hole is irrelevant for the computation of P(BQ; see Figure[2]for an 
illustration. Thus, we can ignore all segments in P * Q that are induced by the 
hole when computing P © Q. It implies that the hole can be removed (that is, 
filled up) before the main computation starts, regardless of the approach that 
one uses to compute the Minkowski sum. 



Figure 2: A small hole H is irrelevant for the computation of P © Q as adding dH and 
7 C Q fills up any potential hole in P © Q related to H. 


Definition 1. A hole H of polygon P leaves a trace in P © Q, if there exists a 
point r = p + q £ d(P © Q), such that p £ dH and q £ dQ. We say that r is a 
trace of H. Conversely, we say that a hole H is irrelevant for the computation 
of P © Q if it does not leave a trace at all. 

Lemma 1. If H leaves a trace in P®Q at a point r, then r is on the boundary 
of a hole H in P © Q. 

Proof. Consider the point r = p+q, which is on the boundary of P®Q, such that 
p £ dH and q £ dQ. Since the polygons are closed, for every neighborhood of r 
there exists a point r' £ P © Q, see Figure [3] Consequently, its corresponding 
point p' = r' — q , which is in the neighborhood of p must be in H. Thus, r' 
must be enclosed by dH © q , implying that r' is inside a hole of P © Q. □ 



Figure 3: Hole H leaves a trace in P (B Q at point r, which must be on the boundary of 
some hole H in P © Q\ see Lemma |T] 

Lemma 2. Let H be a hole in P © Q that contains a point r = p + q £ dH, 
such that p £ dH and q £ dQ; that is, r is a trace of H. Then Vz £ H and 
\/y £ Q, it must hold that z £ H (By. In other words, H C DvyeQ H (By. 
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Proof. As in Lemma [lj there is a point r' in the neighborhood of r, which is 
enclosed by dH © q. Furthermore, there exists a continuous path from r' £ H 
to any z £ H, which means that every z £ H is also enclosed by dH © q, or in 
other words: z £ H © q. 

Now, assume for contradiction that there is a point y 0 £ Q , for which z is 
not in H (B yo- Consider the continuous path 7 that connects q and yo within 
Q. Observe that z £ H © q and z ^ H (B yo are equivalent to z — q £ H and 
z — yo H, respectively. This means that z — yo is either in the unbounded face, 
or in some other hole in P. Now observe that the path z©(— 7 ) connects z—q and 
z — yo■ Thus, since 7 is continuous, there must be a point y' 0 £ 7 C Q, for which 
z — y'o £ P. Hence, z £ P © y’ 0 , which implies z £ P ® Q —a contradiction. □ 

Corollary 1. Let H be a hole in P®Q with r £ dH being a trace of H. Then 
Ms £ dH it holds that s is a trace of H. 

Proof. By Lemma [5] FT C Hv^gQ H ®V □ 

Theorem 1. Let H be a closed hole in polygon P. H is irrelevant for the 
computation of P © Q iff there is a path contained in polygon Q that does not 
fit under any translation in —H. 

Proof. We first show that H is irrelevant for the computation of P © Q if there 
is a path 7 C Q that does not fit under any translation in —H. Assume for 
contradiction that H leaves a trace in P © Q\ that is, there is an r = p + q £ 
d(P © Q), such that p £ dH and q £ dQ. By Lemma [I] the point r is on the 
boundary of a hole H in P © Q. By Lemma [ 2 ] for any point x £ H it must 
hold that x £ H © y My £ Q. Specifically, it must hold My £ 7 C Q. This is 
equivalent to y £ (x®—H) for all y £ 7 , stating that 7 fits into —H under some 
translation—a contradiction. 

Conversely, if there is no path that does not fit into —H then all paths contained 
in Q fit in —H. Thus, also Q itself fits in —H under some translation x with 
x®Q C —H. In this case x + q £ —H for all q £ Q, which is equivalent to —x £ 
H © q for all q £ Q. This implies that —x £ P © Q, whereas — x £ ( P\JH)®Q , 
that is, H is relevant for P © Q. □ 

Corollary 2. If the closed axis-aligned bounding box Bq of Q does not fit under 
any translation in the open axis-aligned bounding box Bh of a hole H in P, then 
H does not have a trace in P © Q. 

Proof. W. 1. o. g. assume that Bq does not fit into Bh with respect to the x- 
direction. Consider the two extreme points on dQ in that direction and connect 
them by a closed path 7 , which obviously does not fit into —H, as it does not 
fit into Bh- □ 

Theorem 2. Let P and Q be two polygons with holes and let P' and Q' be 
their filtered versions, that is, with holes filled up according to Corollary^ with 
P © Q = P' © Q' ■ Then, at least P' or Q' is a simple polygon. 
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Proof. Note that if Bq does not fit in the open axis-aligned bounding box Bp 
of P, it cannot fit in the bounding box of any hole in P, implying that all holes 
of P can be ignored. Since for any two bounding boxes either Bq (jf Bp or 
Bp <£_ Bq holds, we need to consider the holes of at most one polygon. □ 

Consequently, we can remove all holes whose bounding boxes are, in x- or y- 
direction, smaller than, or as large as, the bounding box of the other polygon, as 
an initial phase of all methods. With fewer holes, convex decomposition results 
in fewer pieces. Moreover, when all holes of a polygon become irrelevant, one 
can choose a decomposition method that handles only simple polygons instead 
of a decomposition method that handles polygons with holes, which is typically 
more costly. As for the convolution approach, the intermediate arrangements 
become smaller, speeding up the algorithm. 


3 Implementation 

The software has been developed as part of the 2D Minkowski Sums package of 
CGAL [ZB], and it uses other components of Cgal [Z3]. As such, it is written 
in C++ and rigorously adheres to the generic-programming paradigm the EGC 
paradigms. In the following we provide some details about each one of the new 
implementations. 

3.1 Reduced Convolution 

We compute the reduced convolution set of segments filtering out features that 
cannot possibly contribute to the boundary of the Minkowski sum (see Sec¬ 
tion and in particular complete holes (see Section [2t. Then, we construct 
the arrangement induced by the reduced convolution set j Finally, we traverse 
the arrangement and extract the boundary of the Minkowski sum. We apply 
two different filters to identify valid holes in the Minkowski sum: (i) We ignore 
any face in the arrangement the outer boundary of which forms a cycle that is 
not properly oriented, as suggested in [2]. (ii) We ignore any face /, such that 
(—P ® x) and Q collide, where x £ f is a sampled point inside /, as suggested 
in [l4j. We use axis-aligned bounding box trees to expedite the collision tests. 
After applying these two filters, only segments that constitute the Minkowski 
sum boundary remain. 

3.2 Decomposition 

Vertical decomposition m (a.k.a. trapezoidal decomposition) and triangula¬ 
tion [3] have been extensively used ever since they have been independently 
introduced a long time ago. We provide a brief overview of these two structures 
for completeness and explain how they are used in our implementations. 

5 Currently, we use a single arrangement and do not separate segments that originate from 
different boundary cycles in the summands (exploiting Corollary [I]. We plan to apply this 
enhancement in the near future. 
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Figure 4: Convex decomposition. |(a)| A polygon with holes 
of the polygon in |(a)| |(c)|Triangulation of the polygon in (a) 


(b) Vertical decomposition 


Vertical decomposition for a planar subdivisions is the partition of the (al¬ 
ready subdivided) plane into a finite collection of pseudo trapezoids. Each 
pseudo trapezoid is either a trapezoid that has vertical sides, or a triangle (which 
is a degenerate trapezoid). Given a polygon with holes, we obtain the decom¬ 
position as follows: At every vertex of the polygon, we extend a ray upward 
if it does not escape the polygon, until either another vertex or an edge is hit. 
Similarly, we extend a ray downward; see Figure |4b| In our implementation we 
exploit the vertical decomposition functionality provided by the Cgal package 
2D Arrangements m- 

A Delaunay triangulation for a set of points in a plane is the partition of 
the plane into triangles, such that no point in the input is inside the circumcir- 
cle of any triangle in the triangulation. A constrained Delaunay triangulation 
is a generalization of the Delaunay triangulation that forces certain required 
segments into the triangulation. Given a polygon with holes we obtain the con¬ 
strained Delaunay triangulation confined to the given polygon and provide the 
polygon edges as constraints; see Figure [4c] In our implementation, we use the 
2D Triangulations [28] Cgal package. 


4 Experiments 


We have conducted our experiments on families of randomly generated simple 
and general polygons from AGPLib [4]; examples are depicted in Figure 5a 
and |5b[ respectively. All experiments were run on an Intel Core 2 Duo P9600 
CPU clocked at 2.53 GHz with 4 GB of RAM. For each instance size the dia¬ 
grams in the figures show an average over 10 runs on different input. Every run 
was allowed 20 minutes of CPU time and aborted when it did not finish within 
this limit. 

First, we compared the running time of the implementations of all methods 
for simple polygons available in Cgal (for details, see |8J Section 9.1.2]), the 
new implementations, and Behar and Lien’s implementation; see Figure |7a| 
The reduced convolution method consumed about ten times less time than the 
full convolution method for large instances, whereas the decomposition methods 



















(a) 


(b) 

Figure 5: Randomly generated polygons: (a) simple polygon with 200 vertices, and (b) 
general with 200 vertices and 20 holes. 


were the fastest for instances larger than 150 vertices. 

Secondly, we compared the running time of the implementations of the three 
new methods (i.e., the reduced convolution (RC), the triangular decomposition 
(TD), and the vertical decomposition (VD)) and Behar and Lien’s implemen¬ 
tation on instances of general polygons with n vertices and n/10 holes; see 
Figure [7b} For each pair of polygons, one was scaled down by a factor of 1000, 
to avoid the effect of the hole filter in this experiment. For all executions, 
the reduced convolution method consumed significantly less time than the two 
decomposition methods. Behar and Lien’s implementation generally performs 
worse than our reduced convolution method. 

In order to demonstrate the effect of the hole filter, we compared the run¬ 
ning time of the implementations above fed with a square of varying size (see the 
horizontal axis in Figure[7c]and |7d[ ) and with randomly generated polygons hav¬ 
ing 2000 vertices and 200 holes. Without the hole filter the running time of the 
reduced convolution method increases as the square grows due to an increase 
of the complexity of the intermediate arrangement. Behar and Lien’s imple¬ 
mentation exhibited constant running time, as it performs pairwise intersection 
testing. When applying the hole filter to our methods, the reduced convolution 
method consumed less time than all other methods. The two diagrams clearly 
show the impact of filtering holes. 




Figure 6: Letters from the font Tangerine used for the real-world benchmark, displayed 
with their offset versions, (a) Lowest-resolution “M" with 75 vertices (b) Highest-resolution 
"A” with 8319 vertices. 
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(c) MS of general polygon (200 holes, 
2000 vertices) and growing square ( x- 
axis)-—without hole filter 



(e) MS of a fixed-size circle and an “M” 
with varying vertex count (cc-axis). 




(d) MS of general polygon (200 holes, 
2000 vertices) and growing square (cc- 
axis)—with hole filter 



(f) MS of a fixed-size circle and an “A” 
with varying vertex count (cc-axis). 


Figure 7: Time in seconds for different methods to compute Minkowski sums for two 
polygons, cr-axis: ^vertices of each input polygon, if not stated otherwise. Legend: 
(RC) reduced convolution; (FC) full convolution; (TD) constrained triangulation decomposition; 
(VD) vertical decomposition; (SD) small-side angle-bisector decomposition; (OD) optimal convex 
decomposition; (HD) Hertel-Mehlhorn decomposition; (GD) Greene decomposition; (BL) Behar 
and Lien's reduced convolution. 
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Note that the polygons used for the benchmarks above do not represent a 
real-world case. Instead, the complex shapes rather constitute a worst-case, 
as most segments intersections are inside the Minkowski sum anyway. For a 
more realistic scenario, consider a text, which we want to offset (for example, 
for printing stickers). In Figure 7e we show the running times of the methods 
available for simple polygons when calculating the Minkowski sum of a letter 
“M” (Figu re [6| with a varying amount of vertices and a circle with 128 vertices. 
In Figure [7f| we show the running times of the methods available for general 
polygons when calculating the Minkowski sum of a letter “A” and the same 
circle. For both letters, our implementation of the reduced convolution is at 
least 5 times faster than all other methods. 


5 Conclusion 

All new implementations introduced in this work will be available as part of the 
next public release of CGAL, which now also supports polygons with holes. The 
decomposition approaches that handle only simple polygons outperform the new 
reduced convolution method (which, naturally handles also simple polygons) for 
instances of random simple polygons with more than 150 vertices. However, 
these rather chaotic polygons somewhat constitute the worst case scenario for 
the reduced convolution method. In all other scenarios, the reduced convolution 
method with hole filter outperforms all other methods by a factor of at least 5. 
Consequently, this is the new default method of CGAL to compute Minkowski 
sums for simple polygons as well as polygons with holes. 
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